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We pursue a scenario where the lighter top squark (stop) mass is accessible for the Large Hadron 
Collider (LHC) in the near future, while gluinos and first two generation squarks are too heavy. 
At i/s = 8 TeV, we investigate the identification of stops which decay predominantly into a top 
quark and the stable lightest supersymmetric particle. We use a simple kinematical variable, M3, 
to reconstruct two top quarks which are pair-produced from the stops, in the fully hadronic channel. 
The dominant Standard Model (SM) background for this signal stems from ti plus jets, with one 
top quark decaying into blu, where the lepton is undetected and the v produces missing transverse 
momentum. The lepton identification efficiency is thus crucial in order to estimate the background 
correctly. We identify kinematical variables to reduce the SM background. We find that it is possible 
to achieve signal and background cross-section at similar levels for stop masses around 350 — 500 
GeV for a x° mass of 100 GeV. 



I. INTRODUCTION 

Experiments at the Large Hadron Collider (LHC) have 
recently reported preliminary evidence for a Higgs-like 
particle, with mass in the region of 125 GeV [HI2]. 

In the Standard Model (SM), higher loop corrections to 
the Higgs mass are quadratically divergent; to avoid fine- 
tuning, new physics should appear around a scale of 0(1) 
TeV and cut off the divergences. The problem is most 
severe in the case of the one-loop correction from the 
top sector, since other contributions to the Higgs mass 
are suppressed by gauge or smaller Yukawa couplings. 
Thus, reducing fine-tuning in the SM leads minimally to 
the conclusion that there should a partner for the top 
quark in the sub- TeV regime which is responsible for the 
cancellations. 

The most widely studied mechanism for canceling the 
divergences is supersymmetry, and in particular the dan- 
gerous top quark loops are canceled by the scalar super- 
partner of the top quark, called a stop (t). While the LHC 
has already started to put constraints on the first two 
generation squarks-gluino plane, the paucity of strong 
constraints on stop masses is due to the small production 
cross-section of stop pairs and a huge background from 
top quark production. The current constraints are: (a) 
for a directly produced stop going to top quark plus next 
to lightest neutralino, with neutralino decaying to grav- 
itino plus Z, ATLAS excludes stop masses upto 240 — 330 
GeV using 2.05 fb _1 [3] (b) at the Tevatron, the stop 
mass constraint is about 150 — 180 GeV [HE]. 

In this paper, we will probe a technique for stop 
searches in the following decay mode 

t -> t + X? , (1) 

where the Xi is the lightest neutralino, which we will take 
to be the lightest supersymmetric particle (LSP). In R- 
parity conserving models, the LSP is the main source of 
missing energy in the event. We will always be speaking 



about the lightest top quark superpartner i±, which we 
will hereafter call t. We will not make any assumptions 
about the spectrum, except that the above decay mode 
is kincmatically allowed and dominant. 

The main challenge in such searches is the fact that 
the LHC is a top quark factory and distinguishing top 
quarks produced from stop decay, as opposed top quarks 
produced directly, can be very challenging. There are 
several established techniques of probing the ti system 
or identifying top quarks: 

(i) cut, the Razor, and Mt2'- These techniques rely on 
the identification of two hemispheres to maximize sensi- 
tivity in searches for a pair of heavy colored objects and 
their cascade decays. These are inclusive searches that do 
not rely on the reconstruction of top quark(s); they have 
larger signal acceptance, but larger background [SHTT]. 

(ii) Methods that reconstruct a top quark or two top 
quarks, followed by cuts on topology or kinematics to 
reduce the background. These methods have lower back- 
ground, and lower acceptance of signal. 

We choose to explicitly reconstruct the two-top quark 
system. Within this class of searches, several options are 
available: 

(iia) M3: This is the invariant mass of trijet com- 
binations with highest vectorially summed px- M3 has 
been used in top quark studies at CMS [T2] and CDF 

m . 

(iib) Kinematic Fits: One minimizes y 2 in dijet and 
trijet invariant masses, matching to W and top quark 
masses; for example |15j : 

(iic) Top Taggers: See [16] for a recent review; 

(iid) Bi-Event Subtraction Technique (BEST): This 
method has been used to identify top quark system in 
|14j . Jet combinatorics is reduced by mixing events in 
hadronic decay chains. 

For the two top quark system in Eq. [I] the highest pr 
jets are most likely to be from the top quark, if it is sig- 
nal. Intuitively, M3 should work well in such a system 
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(note that such an assumption is untrue in a SUSY envi- 
ronment with multiple superpartners undergoing decays, 
where the highest pt jets are probably not from the top 
quark; BEST works better there). 

In [XT] [TS] , top taggers have been used to reconstruct 
top quarks coming from stop decay, while in |19j . shape 
analyses of ^-related distributions have been used to 
probe the stop system when the stop-LSP mass difference 
is degenerate with the top quark mass. 

At the LHC, it is expected that the existence of stops 
will be indirectly established initially using inclusive jets 
+ single lepton + fj analysis, however once any excess 
is observed the direct evidence of the stop can be estab- 
lished though the existence of top quarks in the signal. 
Our goal in the paper is to establish the existence of 
two top quarks in the final states along with the miss- 
ing energy in all hadronic channel. In order to make our 
analysis realistic we use PGS4 detector simulation [20] 
and consider W + n jets, Z + n jets and ti + n jets, with 
n < 6 as well as single top + jets backgrounds. Interest- 
ingly we noticed that tt + (3 — 6) jets contribution to the 
background is comparable to ft + (1 — 2) jets. 

Our finding is that simple kinematical selections with 
the M3 variable is an effective tool for stop searches. At 
y's = 8 TeV, we achieve background and signal cross- 
sections at comparable levels for stop masses around 
350 - 500 GeV. 

The outline of this paper is as follows. In Section [TTJ 
we outline our search strategy. 



Schematics of Signal 



Schematics of Background 



III 



we give 



In Section 

our results in detail. The results are also summarized in 
Table |V| appearing at the end of the paper. We end with 
our conclusions. 



II. SEARCH STRATEGY 

We consider the fully hadronic mode 

pp ft* -+ (tx?)(<xS) (bm°x)(bjjXi) (2) 

We investigate samples with at least four non 6-tagged 
jets and at least two 6-tagged jets, along with large miss- 
ing energy. Signal events are generated with ISA JET J5TJ 
+ PYTHIA background events are generated with 

ALPGEN [23] + PYTHIA, followed by PGS4 detector simu- 
lation. 

The main source of missing energy for the tt back- 
ground are neutrinos coming from the leptonic decay of 
a W, while for the signal the dominant source of missing 
energy is the neutralino. Clearly, after the missing energy 
cut, the most critical factor affecting the discrimination 
of signal over background in the fully hadronic mode is 
the lepton veto efficiency. Due to imperfections of the 
lepton veto, ft events with leptonic W decay could be a 
dominant source of background. 

Our method is to (1) reconstruct a top quark using the 
trijet invariant mass M3, (2) use kinematic correlations 
between the constituents of the two (bjj) systems and fx 
to improve the reconstruction of the pair of top quarks 




v j' Lost Lcpto: 



FIG. 1: [Left] schematic diagram of the signal. The stop pair 
gives rise to tt and neutralinos, which are the main source of 
pr- In the fully hadronic mode, the top quarks decay into 
trijet systems. "System A" is the trijet system containing the 
leading pr jet and reconstructed using M3, while the remain- 
ing jets are called "System B". [Right] tt background after 
lepton veto where the lepton is undetected. The main source 
of f)r here is the neutrino from W decaying leptonically. The 
associated lepton passing the veto is termed a "lost lepton" . 



and (3) finally apply M3 again to identify the second top 
quark. We describe these steps below, before showing our 
results in the next section. 

(1) We use M3 twice. First, combinations of three jets 
are made in the sample, keeping one 6-tagged jet and two 
untagged jets in each trijet combination. Next, the trijet 
combination with the largest vectorially summed trans- 
verse momentum Pj^bj^ ^ s chosen. The invariant mass 

It 



of this trijet combination is defined as M3(p l ^, a b d ) 1 ^ 6 ). 
approximates the mass of the hadronically-decaying top 
quark. Similarly, we find a 2 nd leading trijet combina- 
tion Pxbjj an d invariant mass M3(p|,y). Associated 
with M3, we also define M2, which is the invariant mass 
of the two untagged jets in the trijet M3 combination. 

Using A/3, we identify a first top quark, which we call 
"System A" . This is done by calculating \ 2 for the trijet 
and dijet combination corresponding to the leading pt 
combination M3(p£ b j? K ) and also for the 2 nd leading 

combination M3(p|n b • ■), with a mean top quark mass of 
170 GeV and width of 15 GeV, and a mean W mass of 
80 GeV and a width of 10 GeV. The combination with 
the lowest x 2 is then taken to represent System A. We 
call this combination Af3 min . 

We note that this x 2 analysis allows for more signal 
events in the identification of System A. We show the 
results of this analysis in Section |III C[ 

(2) After the identification of System A, we classify the 
remaining b-jet and non-tagged jets to be "System B"; 
thus, we would denote them as (6bJbJ_b). We employ 
various cuts on azimuthal angles between jets and fr. 
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and Mt between bs and pr- These are motivated by 
the fact that for signal, the main source of missing energy 
is the neutralino, while for the tt background, the main 
source of missing energy is the neutrino coming from the 
leptonic decay of the W, or from jet mismeasurement. 
Thus, for example, for the background, pr is aligned 
along bs, as is clear from the schematic diagram shown 
in Figure [T] For the stop decay, however, the correlation 
between the frr i n the form of neutralino and the bs is far 
weaker. The results of this analysis are shown in Section 

mm 

(3) At the final stage, we employ M3 again to iden- 
tify the second top quark, System B. The result of this 
analysis is shown in Section fill E| 



C. M3: Tagging Top System A 

In this section, we use M3 to tag the top quark in Sys- 
tem A, after a pr cut to further reduce SM background. 
The value of the pr cut is determined by maximizing 
the significance for each choice of mass. This is shown in 
Table IE 

TABLE II: pr cut for various choices of masses. All masses 
are in GeV. 



t 

X? 


350 400 
100 100 


450 
100 


500 
100 


550 
100 


400 400 
150 200 


fk Cut 


145 170 


195 


195 


195 


170 100 



III. RESULTS 

In this section, we describe our selection criteria and 
the cross sections after every stage of cuts (see Tabic IV| . 



A. Background 

We generate the following backgrounds with ALPGEN+ 
PYTHIA+ PGS4: tt + n jets with n < 6, single top + 
jets, W (rightarrowrv) +n jets and Z(— > w, tt) +n jets 
with n upto 6. The background cross-sections are given 
in Table [I] While Ref. [T7] considers background with 
ti+ (<)2 jets in the stop analysis, we find that ti + (3 — 6) 
jets contribution is comparable to tt + (1 — 2) jets. 

TABLE I: Main sources of background. "Others" includes 
single top + jets, W + n jets and Z + n jets with 1 < n < 6. 
All cross-sections are in fb. 



Background 


tt + (< 2)j 


tt + (3 - 6)j 


Others 


Cross-section 


2.0 x 10 5 


0.24 x 10 5 


2.8 x 10 6 



B. Baseline Cuts: 6 jets, Lepton veto, pr 

Our baseline selection cuts are: 

(i) A r nonb _j Cts > 4, and at least two loosely tagged 
&-jets. 

(ii) The leading jet has p T > 100 GeV in \r)\ < 2.5, 
and all other jets have pt > 30 GeV in \rj\ < 2.5. 

(m) Lepton veto: We reject isolated electrons and 
muons with pt > 10 GeV in \r)\ < 2.5. The isolation 

criteria are X>tTo < 5 GeV with AR = °- 4 - 

(iv) t veto: We also reject any hadronically decaying 
r with pr > 20 GeV in \rj\ < 2.1. We assume a r identi- 
fication efficiency of 60% and a fake rate of 2%. 

(v) pr > 100 GeV. 



As described in Section [TTJ we identify System A by 
using M3. Figure [2] shows the comparative distributions 
of M3 min and M3(p^^ s ). We improve the top tagging 
by approximately 30% in signal events in the top quark 
mass region. 




M3 (GeV) 

FIG. 2: Distributions of A/3 min and M3(p!^ ns ). The inset 
shows the distribution after M2 mass window cut. A gain of 
~ 30% in signal is obtained by using M3 mm . The luminosity 
is 50 fb -1 . 

We next perform a W mass window cut on M2 min , 
taking 40 GeV < Af2 min < 120 GeV and a top quark 
mass window cut on M3 min , taking 120 GeV < Af3 min < 
220 GeV. 

We show the M3 min distribution after A/2 min mass cut 
in the inset of Figure [2j 

We now proceed to probe the constituents of the "other 
top quark" in System B. 

D. Angular and Mt Cuts: Kinematic Correlations 
between pr and Jets 

We denote the remaining 6-jet and non-tagged jets as 
(^bJbJb)- We clean up the system with various angu- 
lar and Mt cuts, as mentioned in our search strategy in 
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Section [II] The cuts values arechosen based on Figures [3] 
and|U 

(i) A<^(6 B) |? r ) > 1.2 and A^Cjb(i,2),|t) > 0.7, 
where Jb(i,2) refer to the first and second leading jets 
in System B, respectively. 

(m) Mx(bB, Pr) '■ We choose optimal cut values for 
different masses (see Table III I . 



TABLE III: Mt^b , pr) cuts for various choices of masses. 
All masses are in GeV. 



t 


350 400 450 


500 


550 


400 


400 


X? 


100 100 100 


100 


100 


150 


200 


M T {b B ,fk) 


145 155 165 


165 


165 


155 


155 



After the above cuts, we revert to the trijet {bAjAjA) 
in System A with similar angular cuts between missing 
energy and the 6-tagged jet as well as non-tagged jets. 
These angular cuts are efficient in reducing events with 
"lost leptons" . The cuts are chosen based on Figure [5j 

(m) A<t>(b A ,pr) > 1.2 and A0(j A (i,2),Jt) > 0.7, 
where Ja(i,2) refer to the first and second leading jets in 
System A, respectively. 



E. M3 : Tagging Top System B 

As a last step, M3 is applied in System B, followed by a 
W mass window cut on M2 (40 GeV < M2 < 120 GeV). 
The M3 distribution is shown in Figure [6] Our final 
results with 110 GeV < M3 < 230 GeV are tabulated in 
Table EO 

We note that for the point (mj = 350 GeV, m^o = 

100 GeV) we additionally impose A^{b A ,B,fk) < 2.7. 
Also, for the point (m~ t = 400 GeV, m^o = 200 GeV), the 
W mass window cut on M2 of System B was taken as 
60 GeV < M2 < 100 GeV, while the top quark mass 
window was taken as 140 GeV < Af 3 < 200 GeV. 



IV. SUMMARY AND CONCLUSION 

In this paper, we have explored a search strategy for 
a light stop at the LHC at 8 TeV, using M3 variable in 
the fully hadronic channel. The gluino and the first two 
generation squarks are assumed to be too heavy to be 
produced significantly at the LHC. 

We first performed M3 mm to identify a top quark sys- 
tem (System A). Next, we performed M3 again to iden- 
tify the second top quark (System B), along with a se- 
ries of kinematical cuts to reduce the SM backgrounds. 
Throughout this study, we used PGS4 detector simula- 
tion and considered W + n jets, Z + n jets and ti + n 
jets, with n < 6 as well as single top + jets backgrounds. 
Interestingly we noticed that ti+ (3 — 6) jets contribution 
to the background is comparable to it + (1 — 2) jets. 
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FIG. 3: Distributions of A0(6 S , fir), A<p(j m , fir), A0(j s2 , |t) 
for tt background and signal (mj = 400 GeV, m^o = 

lOOGeV). We cut at A(j){b B ,frr) > 1.2, A<p{j B(h2 ), fk) > 
0.7. Here, &s,jsi, and j B 2 denote the 6, leading jet, and 
next leading jet of System B. The luminosity is 50 fb _1 . 



Tables |IV| and [V] are a summary of the search perfor- 
mance for various choices of stop and neutralino masses. 

In summary, a simple kinematical reconstruction tech- 
nique using the M3 variable is an effective tool for stop 
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M T (b B , £ T ) (GeV) 

FIG. 4: Distribution of Mr(bB, fk) for it background and 
signal (m t - = 400GeV,m s o = 100 GeV). We cut at 

M T {b B ,pr) > 155 GeV. The luminosity is 50 fb" 1 . 

TABLE IV: Final significances for various choices of masses. 
All masses are in GeV. The luminosity is 50 fb~ . 



t 


350 


400 


450 


500 


550 


400 400 


X? 


100 


100 


100 


100 


100 


150 200 


s/Vb 


1.29 


1.71 


1.39 


0.81 


0.35 


0.94 0.47 



searches. We find that at y/s — 8 TeV it is possible to 
reduce background down to a level of signal cross-section 
for stop masses around 350 — 500 GeV for a Xi mass of 
100 GeV. 
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TABLE V: Summary of effective cross sections (fb) for stop pair production and the SM background events in our stop search 
feasibility study. Masses and momenta are in GeV. "Other" sources of background include single top + jets, W + n jets and 
Z + n jets with 1 < n < 6. The significance is given at 50 fb ~ . 



Signal tt + n(< 2)jets tt + n(> 3)jets Others 



350 
= 100 



Ang 



Initial 
Baseline Cuts (Sec 

pi > 145 GeV 
System A: M3 (Sec 
ular and Mt cuts (Sec 
System B: M3 (Sec " 
A(j>(b A ,B,pr) < 
Significance (S/ 



IIIBi 



IIIC 

m 

) 



in 

TT 



D) 



760 
9.96 
6.79 
2.65 
0.55 
0.25 
0.14 



2.0 x 10 5 
192 
82.2 
25.8 
1.61 
0.40 
0.24 



0.24 x 10 J 
147 
69.1 
15.6 
1.71 
0.47 
0.25 



B) 



1.29 



0.24 x 10 b 
147 
47.0 
9.73 
1.06 
0.28 



2.8 x 10° 
31.7 
15.8 
3.14 
0.72 
0.20 
0.10 



400 
= 100 



Ana 



Initial 
Baseline Cuts (Sec 

pi > 170 GeV 
System A: M3 (Sec. 
gular and Mt cuts (I 
System B: MS (Sec. 



IIIBi 



Significance (S/VB) 



IIIC 


ec. 


HID 


Ill 





337 
5.55 
3.62 
1.46 
0.44 
0.20 



2.0 x 10 s 
192 
53,4 
15.4 
0.96 
0.26 



1.71 



2.8 x 10 5 
31.7 
11.1 
1.82 
0.54 
0.14 



Initial 
Baseline Cuts (Sec 
450 fk > 195 GeV 

100 System A: M3 (Sec 



IIIBi 



Angr 



System B: M3 (Sec. 



Significance (S/VB) 



IIIC 


cc. 


If ID 


Ill 





160 
2.52 
1.61 
0.62 
0.25 
0.12 



2.0 x 10 s 
192 
34.5 
9.17 
0.55 
0.17 



0.24 x 10 J 
147 
31.9 
6.32 
0.69 
0.14 



1.39 



rr 



2.8 x 10 6 
31.7 
8.08 
1.30 
0.40 
0.06 



2.8 x 10° 
31.7 
8.08 
1.30 
0.40 
0.06 



Initial 
Baseline Cuts (Sec 
500 fk > 195 GeV 

100 System A M3 (Sec 



IIIBi 



Angv 



System B: M3 (Sec. 



Significance (S/VB) 



IIIC 


Jl 


Sec. 


[11 D 


1111 


3p 



80.5 
1.21 
0.86 
0.32 
0.15 
0.07 



2.0 x 10 
192 
34.5 
9.17 
0.55 
0.17 



0.24 x 10° 
147 
31.9 
6.32 
0.69 
0.14 



0.81 



3~ 



2.8 x 10° 
31.7 
8.08 
1.30 
0.40 
0.06 



m-o 



Initial 
Baseline Cuts (Sec 
550 pi > 195 GeV 

100 System A: M3 (Sec 



IIIBi 



Angr 



System B: M3 (Sec. 



Significance (S/VB) 



IIIC 


cc. 


If ID 


Ill 





43.0 
0.57 
0.43 
0.14 
0.07 
0.03 



2.0 x 10 
192 
34.5 
9.17 
0.55 
0.17 



0.24 x 10 J 
147 
31.9 
6.32 
0.69 
0.14 



0.35 



Initial 
Baseline Cuts (Sec 
400 pi > 170 GeV 

150 System A: M3 (Sec 



IIIBi 



Angr 



System B: M3 (Sec. 



Significance (S/yB) 



III 


1 


Sec. 


HID 


Ill 





337 


2.0 x 10 b 


0.24 x 10 b 


2.8 x 10 b 


4.78 


192 


147 


31.7 


2.76 


53.4 


47.0 


11.1 


1.01 


15.4 


9.73 


1.82 


0.23 


0.96 


1.06 


0.54 


0.11 


0.26 


0.28 


0.14 






0.94 




337 


2.0 x 10 b 
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FIG. 5: Distributions of A</>(b A ,fk), A<p(j A i,fk), and 
A(/>{jA2,Pr) for tt background and signal (m t - = 
400GeV,m s o = 100 GeV). We cut at A0(6 A ,|r) > 1.2 and 

A0(iA(i,2), Pr) > 0.7. Here, b A ,jAi, and j A 2 denote the b, 
leading jet, and next leading jet of System A. The luminosity 
is 50 fb -1 . 
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FIG. 6: Distribution of M3 in System B, after requiring 
40 GeV < M2 < 120 GeV. Displayed are: "other" sources 
of background (single top + jets, W + n jets and Z + n 
jets with n < 6), total background including tt + n jets, 
and total background plus signal for our reference point 
(m~ = 400GeV,m s o = 100 GeV). The luminosity is 50 fb -1 . 
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